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' The nonUnear a model (NLcrM) epitomises a field-theoretical approach to (interacting) electrons 

\ in disordered media. These lectures are aimed at the audience who might have vaguely heard about 

■ the existence of the NLcrM but know very little of what is that, even less so of why it should be used 

04 ' and next to nothing of how it can be applied. These what, why and mainly how are the subject of 

the present lectures. In the first part, after a short description of why to be bothered, the NLcrM is 
derived from scratch in a relatively simple (but still rather mathematical) way for non-interacting 
electrons in the presence of disorder, and some illustration of its perturbative usage is given. In the 
OO ' second part it is generalised, not without some leap of faith, to include the Coulomb repulsion and 

superconducting pairing. 

^ , I. INTRODUCTION 

' Starting from the seminal paper by Wegneri, a field-theoretical description based on the nonlinear a model 
1) ' (NLcrM) has become one of the main analytical approaches to quantum transport and thermodynamics in 
disordered electronic systems [2-24] . Quantum effects are due to interference of electron scattering from different 
^ ; ' impurities. These effects are relevant when the scattering mean free path ^ is small compared to the dephasing 
, length which characterises the scale above which interference is destroyed by an irreversible phase-breaking 
' in electron's wave function due to inelastic interactions. As diverges when T— >0, the quantum interference 
I , effects are always important at low temperatures. Although the quantum corrections to transport coefficients 
■ O ' are small in metals as a power of the small ratio Ap/-^ (Ap is the Fermi wavelength), they are essential in 
electron transport as they govern a nontrivial dependence on temperature, frequency or external field3^^*2&. 
Furthermore, the electron-electron interaction in the presence of the scattering from impurities also results 
in similar (small in Xf/£) corrections both to the transport coefficients and to thermodynamic quantities (see 
Refs.^'''^® for reviews), in a stark contrast to the case of pure normal metals where the interaction manifests 
itself only via the Fermi liquid parameters, which renormalise electrons spectral coefficients. The NLcrM gives 
analytical tools for describing both the weak disorder regime (Ap <C i) and a regime when the disorder increases 
with Xf/£ approaching a critical value of order 1 where electronic states at the Fermi energy are localised^^iSS. 
] There are two main characteristic features in the cr-model approach to the disordered system: (i) one performs 
. the ensemble averaging over all configurations of disorder at the first stage of calculations; (ii) one separates 
' "fast" and "slow" degrees of freedom and deals only with the latter. In the present context, "fast" implies 
m : energies of order ep and momenta of order pp while "slow" means energies smaller than h/rd = hvp/i and 
momenta smaller than h/ Such a scale separation works for a weak disorder. Quite often, the weak disorder 
condition \-p <^ i is formulated as go S> 1 where g is the dimensionless conductance (loosely, the conductance 
measured in units of e^/fi) and cx {e-pTc\Y^^ is the dimensionless conductance of a cube of size l"^ in a 
H _ d-dimensional system. The slow modes in a system of size L ^ £ describe diffusive propagation of electrons, 

. and the standard NLcrM is the theory of interacting diffusive modes (diffusons) 
^ ' The functional of the cr model has the following form: 
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where D = VpTc\/d is the diffusion coefficient, i/q is the one-electron density of states, d is the spatial dimen- 
sionality and Lo is the frequency. The Hermitian matrix field Q{r) (whose symmetry and dimensions depend 
on this or that particular way of the ensemble averaging) obey the following constraints which make the model 
nonlinear: 

Q^ = l, TrQ = 0. (2) 

Writing down the above functional hardly explains for non-initiated what is the meaning of the model. It allows 
me, however, to formulate the main goals of these lectures: 

1. To outline a pedestrian but consistent step- by-step derivation of the above functional for non-interacting 
electrons in a random potential. 
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2. To generalise the model, with the help of one or two leaps of faith, for including the Coulomb repulsion 
and pairing attraction of electrons. 

3. To illustrate using one or two relatively simple examples how the NLcrM works, both for noninteracting 
and interacting electrons. 

But before attending to this programme, it is useful to give some ideas of why the NLaMis to be used at all, or 
where there are its advantages as compared to the usage of other analytical means, like a straightforward dia- 
grammatic approach or a semi-phenomenological random matrix theory (whose applicability to the description 
of non-interacting electrons in disordered systems in a so-called ergodic regime has been rigourously established 
with the help of the NLctM^-*'). 

The first application of the NLcrM was a derivationSiSi^ of the renormalization-group (RG) equations of 
the one-parameter scaling theory-- of the Anderson localisation, illustrated on Fig. ^ The (3 function, (3{g) = 
din g / dlii(L / £) , depends only on the dimensionless conductance g and describes how it changes with scale or 
disorder. For d = 3, conductance increases with the system size for P{g) > and decreases for P{g) < 0, 
becoming an ideal metal or ideal insulator, respectively, in the limit L oo. In this case (3{g) = corresponds 
to the critical point for the metal- insulator transition. For d — 2 and d = 1, P{g) is always negative so that 
no truly metallic state can exist in the thermodynamic limit. However, scaling considerations alone could only 
prove^^ that (3{g) = {d — 2)g for g ^ 1, thus leaving some uncertainty in the marginal d — 2 case. The 
uncertainty has almost immediately been resolved^^ with the help of the diagrammatic perturbation theory by 
proving that g/go = 1 ^ 9o^ InL/i for d = 2, i.e. that the correction to the classical behaviour, P{g) — 0, was 
indeed negative. Assuming also the renormalisibility, supported by the cancellation of all the corrections to g of 
order [g^^ InL/i]'^ also found in Ref.^®, would support the large-g limit of the one-parameter picture on Fig.^ 
which, however, has been confirmed directly only within the NLctM. The result /3{g) = —l/g + o{l/g'^) found in 
[2-5] means that all the main logarithmic contributions, [g^^ \nL/£]"', should be cancelled in the perturbation 
theory, while diagrammatically such a cancellation has only been demonstrated^^ for n — 2. 

Having been originally designed for describing the electron motion in the diffusive regime, the NLcrM has 
later been generalised both for a treatment of the (non-perturbative) ergodic regimei^, where it was used to 
justify the applicability of the random matrix theory for the description of electrons in small metallic grains, as 
had been suggested earlier'^^, and for a new approach to the ballistic chaotic regim9i2ii^iiSii&. Thus, the NLcrM 
covers all the variety of conditions in electronic disordered systems (see Fig.|2l. Moreover, this an indispensable 
analytical tool to describe changes in properties of the disordered system approaching, with decreasing g, the 
region of the metal-insulator transition in d > 2 or of the crossover to the strong localisation in d = 2. 

Further extensions of the NLcrM have been related to the emergence of mesoscopicsiS^i^^iSSiS. When it 
became clear that fluctuations of practically all observable quantities in sub-micron and nano-samples are 
essential so that the ensemble-averaged quantities alone are not sufficient for a full characterisation of disordered 
systems, the NLcrM turned out to be a natural tool for finding entire distributions of different observable 
quantitiesii*i2iiii4, and describing their change with increasing disorder^^. 

The NLcrM has also been extended to include interactionaiSiSiifliiLiSiiSiSfliSi which allowed to reproduce 
earlier perturbative resultsSLSS with a proper account for the Landau Fermi-liquid parameters, and to derive 
the renormalization group equations describing (at least at a qualitative level) a metal-insulator transition 




FIG. 1: Schematic shape of the /3 function in the one-parameter scaling theory 



3 



in disordered interacting systems. The NLcrM has been apphed also to superconducting systemsiSiSSiSi, in 
particular to include effects of Coulomb interaction, e.g., lowering the transition temperature by disorder - 
in a (well understood) deviation from the Anderson theorem2£; and also to describe inhomogeneous or granular 
superconducting systems. The interest in this approach has been greatly enhanced by the recent discovery^S 
of an apparent metal-insulator transition in 2D disordered systems in zero magnetic field. 

The starting point for all further considerations will be the TOE model given below. First I introduce a 
functional integral representation for a simplest variant of the model that does not include interactions; then 
I will show how to derive the NLctM from such a representation by excluding "fast" degrees of freedom; after 
that I will extend the NLcrM by including the Coulomb and pairing interactions. 



II. FROM THE TOE MODEL TO A FUNCTIONAL INTEGRAL 



Here TOE stands for "Theory of Everything" . In contrast to the high energy physics, there is no problem 
to write down the TOE Hamiltonian in context of electrons in solids: 

i i ij 

Here the second term represents schematically both lattice and disorder potential, while the last represents all 
possible two-particle interactions (for example, the Coulomb repulsion or the Cooper attraction). In order to be 
really the TOE model, the above Hamiltonian should also include spin terms, as later in this paper in relation to 
superconductivity. Naturally, either with or without spin terms, there is no hope for a rigourous approach to, let 
alone the exact solution of, this model for an arbitrary disorder/interaction strength. Moreover, using this model 
for many problems in strongly correlated systems is no more fruitful than starting entirely ab initio. However, 
when the aim is to describe low-temperature properties of interacting electrons in disorder environment, which 
are governed by low-energy modes only, the TOE model is perfectly adequate and one does not need to go 
to a next hierarchial level. This docs not make the model easily treatable, let alone straightforward. It 
became clear about 20 years ago that the interplay of disorder and interactions makes disordered conductors 
drastically different from pure metals, where the interactions reveal itself only in renormalisation of Fermi liquid 
parameters. In the perturbative regime the electron-electron interaction leads to nontrivial corrections to both 
transport coefficients and thermodynamic properties. These corrections govern numerous physical effects, like 
zero-bias anomaly or negative anomalous magnetoresistance, which have been summarised, e.g., in reviewsSLS^. 
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FIG. 2: Different regimes for disordered systems: the ergodic regime corresponds to the times longer than L^/D, the 
time of diffusion through the sample, while the ballistic to the times shorter than Tei; at the critical disorder, <? 1, the 
former regime shrinks to the energies of order or smaller than the level spacing S and the later to the energies of order 
of ef, while the diffusive regimes evolves to a critical one embracing all the energies between S and ep--. 
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Note than in the present context a small perturbation parameter is g <C 1, independently of a value of the 
interaction parameter /Hvy'- in contrast to pure metals the condition of high density (Ap ob = /me^) is 
not necessary for a diffusive weakly disordered system, when Ap, ob ^ ^ <C i. A set of interesting questions is 
related to behaviour of the system when the disorder is increasing. 

The most consistent treatment of the TOE model which includes describing its properties with increasing 
disorder can be achieved within the NLcrM that can be derived from Hamiltonian Q for any reasonable interac- 
tion Vij, including both Coulomb and pairing interactions. However, disorder makes the model nontrivial even 
in the absence of any interaction. So the first task here is a derivation of the NLcM in the most straightforward 
case oiVij=Q. A non-trivial step in this derivation would be the ensemble averaging over all realisations of 
disorder. But before coming to that, it is necessary to represent Green's functions as functional integral. 



A. Introducing Green's functions 



Any physical quantity (any response function) can be represented, at least in the absence of interaction, as a 
product of one-particle Green's functions. To this end, let's first write down the standard formal representation 
for the one-particle Green's function G(e), defined by (e — H)G{e) — 1. For Vij — 0, one formally diagonalize 
the Hamiltonian as H \ a) = |a) . Using the completeness condition, 1"^) ("^1 = ^las (with h = 1 

from now on) 

G^={e-H±iS] ={e-H±iS 



e — Sa^iS 

where infinitesimal S is introduced to avoid divergences at e = Sa- On recognises that are the standard 
retarded and advanced Green's functions^'' which can be re-written in the real space representation as 

GHry> . irlG- K) ^ j: <^^^<'^;^^ . j: <^'r'_^f' ' (5) 

a ^ a ^ 

In the presence of disorder, one does not know either eigenfunctions, (pa{r), or eigenvalues, of H. The task 
is to represent the average products of G without any reference to unknown parameters. To this end, I shall 
first represent G as a "functional integral" which is just a multiple product of Gaussian integrals. 



B. From Gaussian to functional integrals 



The first step is representing the denominator of Green's function as a Gaussian integral. Noting that (with 



In = / da;e*(^ 



i{e+-ea) 

one represents (e^— £q)~^ as {i/tt)!'^. (One assumes from now on that energy is dimensionless, which will make 
all the "action" functionals below also dimensionless). Similar representation would be, of course, possible for 
(e~— Eq)"^ but not just for (e— e^)"^: the convergence in /o is ensured by the presence of iS. Representing Jg 
as a product of integrals over dx and dy and changing variables to c — x + iy and c* = x — iy, one finds 

^^*^%-*(^"--)=^Zo. (6) 



e+ — J 27r 

This gives the following representation of det(e+ — H)^^ = Jlal^^ — £a)~^ as a product of Gaussian integrals: 

where 'Dc*'Dc is a symbolic notation for the measure of integration, which is just a product over all dc*dcQ,. 
Of course, the above determinant is defined as a product of an infinite number of (well-defined) integrals, and 
as such, it may appear to be ill-defined. However, if one does not pursue mathematical rigour, this is quite a 
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convenient representation which will be used as a building block for that for Green's functions, the latter will 
be free from trivial divergences that may appear in calculating det(e — H). 

Let us represent (exactly) the sum in the exponent in Eq. by an integral. Using the orthonormality 
condition for the eigenfunctions of H, one obtains 

^c;(e+ -e„)c„ = ^4(6+ -ea)c^ / d'^ripl{r)(pp{r) 



As c and c* were variables of integration which took all possible values in the complex plane, the complex 
conjugate fields ip and ■(/;* take all the possible values for any argument r. Using the identity (|SJ), the determinant 
of Eq. O can be represented as the following functional integral over the fields ip* and i(j 

Zi = dct(e+-H^ ^ = jvi^*V^e'^'^ , S+ = jd'^rij*{r) [e+ - Hj-<p{r) (9) 

(An index in Zi is introduced here to make some notations in what follows more consistent). The measure 
of integration is written here symbolically as 'D'ip*'D'ip in understanding that, using Eq. JSJ, one can always 
reduce the i/)-integral to a (relatively) well defined integral over c and c*. To illustrate how it works I shall 
substitute this integration measure and the definition of fields ?A and ip*, Eq. ((HJ, into the functional average 
—i {'4'*'4')+ = f e'^^'^'D'il!*'D^p, thus demonstrating that this functional average provides a representation 
for the Green function of Eq. Q : 

dc%dcr. 



-0 unless a— a' 



Zi ^ (e+ - e„)^ -f/ e+ - sp ^ e+ - £q £+ - ep 

a \ "/ B^a ^ a all fl ^ 



G+(r,r'). 



However, the functional integral would be useless, had one every time needed to return to the c representation 
in order to make a calculation. It is sufficient to realise once and forever that the functional integral is a short- 
hand notation for a product of Gaussian integrals, so that all the standard operations for Gaussian integrals 
can be directly translated into the language of functional integrals. To illustrate such a translation, let us note 
that in the above calculation the following elementary integral has been used: 

_ [ dc^dc^ ^ _ i _ i 

where Zq is the integral ©. There are many ways to get such a straightforward result. Let us choose that 
which is directly generalisable to the functional integral, representing Iq above using parametrical derivatives 
with respect to the "source fields" h and h* 

J_ r dc*dc * .c'Ac L_^!_ f '^^*'^^ r(c'Ac+h'c+c'h) 

° Zq J 2m Zo dh* dh J 2m 

1 [ ^{c-Ac-h■'A-^h) ^ ^h'A-^h ^ i_ r,^ 

Za dh* dh J 2m dh* dh A ' ^ ' 

where a shift of variables, c ^ c — A^^h and c* ^ c* — h*A^^, has been used. It was implied that both h and 
h* are put to zero after calculating the derivatives above. Now I apply practically the same procedure to the 
functional integral defined above in order to represent Green's functions. 
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C. Green's function as a functional integral 

To show directly that —i {'ip*'ip)_y. represents Green's function, 1 express it as a functional derivative: 

i 5^ 



-i{r{rimr2))_ 
S+{h,h* 



Zi 5h*{ri)5h{r2) 

ip*{r) (e+ - H^j ij{r) + h*{r)iP{r) + il)*{r)h{r) 



(11) 



Noting that, by definition, \ e'^ — H j = Gj and making a shift of variables similar to that above, tp—G-h 

and ip* ip* — h* -G (with {G-h){r) being a convolution in the r space) one reduces S+{h, h*) in the functional 
integral to the following expression: 



S+{h,h*) ^ S+- j[h*{r)Gt{r,r')h{r')]d'^rd'^r' , 

so that performing the functional differentiation one finds, indeed, that 

1 S^Zi{h*,h) 



G?(n,r.) . JvrVi. rir.mr.)e-^ Sk^in) Shir.) 



Zi{h*,h) = V'lP*V^ Q^Sih'■M) 



h =0 



(12) 



That completes the representation of Green's function as a functional integral. If one thinks that the above 
derivation was not "pedestrian" enough, it would be sufficient to compare it to the standard way of introducing 
functional integrals (see, e.g., Ref4i). Note that (|ll|l is not an action as usually understood: in the absence of 
time integral, it does not even have the dimensionality of /i, and is simply assumed dimensionless. Unfortunately, 
such a simplified way is only possible for a quadratic Hamiltonian which can be straightforwardly diagonalized. 
Only the presence of disorder makes it nontrivial: H in Eq. contains a random potential V . The ensemble 
averaging over the distribution of V would result in an effective functional containing higher powers of tp* and 
-0. If _ff is nonlinear to begin with, like the Hamiltonian of the TOE model Q, one should choose a hard 
way^^ to get a functional representation for Green's functions. However, there is an alternative. After having 
learned a bit how to deal with the functional integrals that emerged for the noninteracting problem, one would 
be able to use some leap of faith to extend considerations to include interactions. Such an approach would not 
be self-contained but would remain, in my opinion, reasonably consistent and believable. 

A certain extension related to future considerations is useful right now. Particle statistics are totally irrelevant 
for one-particle problems, but in extending considerations to the interacting particles one shall be mainly 
interested in fermions. The functional integration described above could be extended to the case of interacting 
bosons rather than fermions as the commuting classical fields in the functional integrals correspond to operators 
with certain commutation relations in the Hamiltonian representation. As fermions are described by operators 
with certain anticommutation relations, it is not that surprising that the functional integral representation for 
fermions must be formulated in terms of anticommuting classical fields. There exists a very solid mathematical 
foundation of the theory of anticommuting classical variables*-- reduced to a physical level of rigour in different 
derivations*-^-^^ of the NLcrM. However, all one needs for the purpose of present considerations is a relatively 
moderate leap of faith that allows one to introduce a fermion representation at a one-particle level. 



D. Fermions do it like bosons 



but upside down: instead of having dot ^ in the partition function for bosons 10, one has det in the 
appropriate fermionic representation: 



Zi+ = det ( £+ - ) = IVipViPe'^* , S+ = 



d'^r 



zt 



VipVil; i/'(ri)0(r2)e*' 



(13) 
(14) 
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where ip and ip are classical anticommuting mutually conjugate fields. The conjugation requires some care. For 
the bosonic classical fields it was implied that ijj,'4'* = o ± ib, with a and b being real fields. For the fermionic 
fields it is more convenient to introduce the conjugation operation in matrix notations: 

"^^Tll^)' ^"72^"''* (^*)* = -^' v^^(c^)^ (15) 

The last equation above defines the so called "charge-conjugation" matrix C (which has nothing to do with 
the electric charge). The symmetry imposed by such a conjugation will not be used but in section Hvl In 
what follows, it is sufhcient to know that the classical fermionic fields anticommute (which implies not only 
>f>f* — ~>t* >c but also >!:^ = >f*'^ = 0) and the partition function (|13|l is given by det(e — ¥1) rather than by 
det"\e-7?). 

Now everything is ready for the ensemble averaging over disorder and the only question to ask is what 
to average. The averaged one-particle Green's function (G) is not an interesting object: in the absence of 
interactions disorder does not effect physical quantities which are directly expressible via (G) , such as density 
of states. Most of response functions of interest such as density-density correlations, conductance, mesoscopic 
fluctuations of any observable, etc., can be normally expressed via the averaged products of retarded and 
advanced Green's functions. So the simplest product to average over disorder is G'^G~_^, with G+ represented 
by Eq. ifT^ and G~ represented as 



G;_^(ri, ra) = -— jViPV^, V(r-i)V(r-2) e^*" , (16) 

where differs from S+ in Eq. H13|) by substituting e — uj^ = e — {lo -\- i5) for = e + i5. In general, one 
averages a closed loop, i.e. a product like G(ri, r2)G(r2, r^) . . . G(r„, ri). This allows one to introduce quite an 
economical representation, using one source field for a product %p(r)il}{r) of two fields taken from the adjacent 
Green's function: 



K^{r,r') = G\ ^,Jr,r')G- ^ .Jr' ,r) ^ -^^^^^^^^ 



(17) 



ri=0 

iSi(ri,iAj) 



Zi{r],uj) = JvipVip e'' 
Siiri,u;)^ y'|v^(r")(^-|-V+ic.+A^tA(r") + (V^r,tA),„|dV'. (18) 



Here ^ = p /2m — ep; the frequency arguments have been shifted for convenience; — ,■)/' j, with ^ 

having been used in representations for G^, respectively; "0 is a corresponding vector column; r\ is an arbitrary 
4x4 matrix field (whose diagonal elements, rf^^ and 77 , are redundant and can be put to zero from the very 
beginning, if convenient, while and r\ ^ arc 2x2 matrices by themselves, as each i\r^ is a two-component 
field, Eq. (|15|l '): the constant matrix A = diag(l, 1, — 1, — 1) is necessary as lo^ enters with the opposite sign 
in the actions Z\ is the functional integral in the r.h.s. of Eq. H17|) calculated at r; = 0. As the 77-part of 
the action is not affected by the averaging, the resulting functional would be applicable to any product of G"*" 
and G~ subject to a proper extension of the matrix space (each G would require its own pair of fields in the 
functional representation) and an appropriate choice of the source fields 77. 

III. AVERAGING OVER DISORDER 

For simplicity, one assumes the random potential V"(r) in Eqs. © and (|18|l to be a Gaussian white noise with 
(F(r)) = and 

^V{r)V{T')) = -^8(r-r'), (19) 



which means that Viy) is taken from the following distribution: 



F{V^(r)} =7Vcxp 



(20) 
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where is a normalisation factor. Although the choice of a white-noise Gaussian potential makes calcula- 
tions considerably easier, a generalisation to a slowly varying potential would not basically change results-'^^-'^ 
provided that the mean free path € is substituted by the transport length £tr- With this choice, the disorder 
averaging (. . .) reduces to taking a functional average with the distribution H2U|) . Had one to average Zi, this 
would be a straightforward Gaussian integral again. It is not straightforward to average InZi, though. If 
one first takes derivatives with respect to the source fields in Eq. H17(l . it is the presence of the denominator, 
17.^1(77 = 0), that makes the averaging tricky. 

This denominator is the partition function given by Eq. @ in the bosonic and by Eq. in the ferniionic 
representation. Comparing these expressions, one can see that using a mixed, supersymmetric fermion-boson 
representation would cancel out this denominator. This idea— turned out to be extremely fruitful and led 
to enormous progress in theory of non-interacting particles in disordered and chaotic systems. Numerous 
achievements of the supersymmetry method have been summarised in a book^ and a recent reviewi^. However, 
there is no simple way of applying this method to a many-particle system of interacting electrons, so that some 
other method is required. 

After a few earlier attempta^^*^^*', it has recently been dcmonstratedi2*iS that the Keldysh technique^S*^ 
provides a viable route to a a field-theoretical description. In this technique Z = 1 by construction, which makes 
the ensemble-averaging straightforward. However, this technique remains quite cumbersome and, although it 
might be indispensable in describing non-equilibrium problems, a considerably simpler alternative exists for 
interacting electrons in disordered systems at equilibrium (i.e. in the linear response regime for all transport 
problems). 

This alternative is a so called replica trick, which is possibly the oldest method ^^i^" used for the disorder 
averaging. It is based on the following identity: 

2^ " — 1 

InZi = lim — . (21) 

To apply it one replicates the fields in Eq. H17f) . 1^ = ("01^) • ■ • i i^n ^i^i > ■ • ■ > V'n )' thus replacing it by the 
following: 

'^"(■■•■■■)^iL".i''{s|f?foL' ^-^.. = ^i' = >™»e-<'). (22) 

The functional integration in Eq. 122|) implies a product of n identical integrals; S'„ is given by the same 
expression (|18|) as 5*1 but in terms of the replicated fields and with A = diag(]l2n, — l2n) being a block-diagonal 
An X 4n matrix; is a block-off-diagonal 4n x 4n matrix with 77++ = r) = 0„ and rj"' = rj ^ being arbitrary 
2n X 2n matrices: 



-i ~V{r) + ^uj+ A + T] 



*(r)U'*r. (23) 



An extra factor of l/2n in Eq. H22|l compared with Eq. (|21|l is due to the fact that the trace there gives the sum 
of 2n identical terms. From now on I will be dealing only with the replicated Z and S and thus omit the index 
n there. 

As the representation H22() does not contain a logarithm, one only needs to average Z over the disorder 
distribution of Eq. H20|l . which, for any integer n, is a straightforward Gaussian integration that involves only 
the V^-dependent part of Z: 





f {■^■^fd'^r 


47ri'oTci J 





(24) 



J\f ^ 

Thus one arrives at the following effective action: 

8(7]) = y"|^(r)(^-e+ic^+A + r7^*(r) + ^-^(**)'|d'^r. (25) 

After calculating K,^{r,r') for any positive integer n one needs to take the replica limit n ^ which is by no 
means a well defined mathematical operation. 

However, as long as one is interested in perturbative corrections only (including the RG corrections), one 
can talk about the replica method, not a trick, as illustrated in Fig. O All the perturbative corrections in 
the standard diagrammatic techniques are just crosses (corresponding to scattering from impurities) connected 
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FIG. 3: Diagrammatic corrections for one-particle Green's function: dashed lines correspond to pair correlations of the 
random potential, Eq. 1191 , or equivalently to the "interaction" (the bilinear term) in Eq. 1251 ; wavy lines correspond to 
a real interaction introduced in section IVTI 




by solid lines corresponding to Green's function, like diagrams (a) - (c), where the dashed lines connecting 
the crosses represent the pair correlations of the random potential (|19|l resulting from the ensemble averaging. 
The field theory with the functional (|25|l with any integer n exactly reproduces all such diagrams but also 
yields those like (d) and (e) with vacuum loops. However, each loop contains a sum of n identical replicated 
contributions which lead to an extra factor of n thus making such diagrams to vanish in the replica limit 
n —^ 0. Note that in the supersymmetric technique such diagrams vanish as each loop is the sum of fermionic 
and bosonic contributions having the opposite signs. However, in the latter approach it would be difficult to 
include diagrams like (f) that contains a loop corresponding to the electron-hole excitation in the presence of 
the interaction, whereas the replica approach would encompass such contributions rather naturally (section 
ED). 

The functional of Eq. H25|l is fully equivalent to the original model of Eq. ||2Jl in the absence of interactions. 
However, there is little to gain from it without further approximations. The most important one is based on a 
reduced description that excludes all details of electronic motion at scales smaller than the mean free path £. 



IV. THE NLcrM FUNCTIONAL 

A. Hubbard— Stratonovich decoupling 

This is the name of the operation performed in order to exclude all "fast" modes of electron motion. Its first 
step is nothing more than the Gaussian integration like that in Eqs. H10() and H24|) . but the integration which 
is taken "backwards". In Eq. (|24|l . one integrates over the random potential V thus arriving at the action (|25|l 
that does not depend on V (and thus contains no randomness) but is no longer quadratic (i.e. integrable) in 
the fermionic fields. In order to decouple the quartic term in Eq. (|25|l . one introduces a new field, Q, which 
allows one to "unwind" the integration in Eq. H24() and represent the quartic term as a functional integral over 
Q, thus making the action quadratic in the fermionic fields again and finally integrating over to arrive at 
the action in terms of the new field Q. Of course, one could perform the ^'-integration from the very beginning 
obtaining the action in terms of the disorder potential V; however, such an action would not be very useful. 
The introduction of the field Q is useful if it exploits fully the configurational space defined by the conjugate 
fields 4* and and excludes fast degree of freedom; to this end, Q is chosen as a matrix field having the same 
rank and symmetry as ^' (X" ^' (the latter is just a notation for an outer product of two vectors that results in 

a, An X An matrix). Noting that tr (g) = — ('I'^) and choosing Q to be dimensionless, one can write the 
following identity that "unwinds" the integration l|24|): 



exp 



1 



A-KVqTcX 



■VQ 



■ exp 



trQ2 



2t, 



dS 



(26) 



The normalisation integral Mq is just the same Q-integral at 'J = = (and it is anyway irrelevant in the 
replica limit) and the Gaussian integration both here and in Eq. (|24|l is performed with the help of a shift of 
variables fully analogous to those in Eqs. H1(J|I and Hll|l . After this, the ^'-dependent part of the action is, indeed, 
quadratic, with the term i^^^Q"^ being added into the brackets in Eq. (|25|) . The Gaussian integration over 
the fermionic field in Eq. H26|) is performed exactly like in Eq. (|13|l (the rank of matrices does not matter) and 
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a 



FIG. 4: Diffuson (a) and cooperon (b) ladders: in diagrammatic techniques either pictures with triple dashed lines or 
those with wavy lines can be used to denote the ladders; the latter are more convenient in the present context as they 
are in direct correspondence with diagrammatic notations in the NLctM, section |3 they should not be confused with 
interaction lines, Fig. 

results, after using the remarkable matrix identity Indet A = trln^, in the effective action which depends 
only on Q : 



UJ - I 

-A + r]-^+—Q 
2 2tc\ 



(27) 



where Jo = ^/{vqL'^) is the mean level spacing. The action has been renamed into T in order to stress that the 
functional averaging in Eq. (|25|l is performed with the Euclidian weight e~'^ rather than with the oscillating 
one, e*"^. The partition function Z, Eq. ^22\ . is now given by 

Z = JvQexp{-J^[Q]) , (28) 

and the measure of integration will be specified later. The matrix Q is Hermitian as it has the same symmetry 
as the product 5* (g) ^ which is Hermitian due to the definition of the fields in Eq. H15|l . This is not the only 
symmetry of Q, as the charge conjugation introduced in Eq. H15|) is essential and explains the appearance of 
the factor 1/2 before the trln in Eq. (|27|) . mysterious for generations of students. 

If a diligent student had taken a course of action described between Eqs. (|26|) and H27(l . he would have duly 
arrived at Eq. (|27l) but without this 1/2. Such a result, moreover, would be formally correct but with a matrix 
Q ill suited for describing the slow variables. The point is that the field Q should be chosen in such a way that 
it would correspond only to the slow part of tr (ji> (g) ^) in the l.h.s. of Eq. (|26|l while the fast Fermi oscillations 
contained in each field -0 were suppressed. Such a slow part could be chosen in two ways: it corresponds either 
to the 5'^' pairing - with a small transferred momentum in the Fourier space, or to the '5 ^E* (and ^'5') pairing 
- with the small sum momentum. The charge conjugation introduced in Eq. (|15|l allows automatically for the 
existence of the two pairing channels, provided that the field Q is chosen to have the same charge-conjugation 
symmetry as the product i8) ^E", i.e. Q = CQ*C^^. Allowing also for the fact that Q is Hermitian (which is, 
of course, necessary for at least a superficial convergency of the Q integral in Eq. H2t)|) ). one can represent the 
An X 4n matrix field Q as 



^^[c+ <F j ' d+=d, = -c, (29) 

where d and c are 2n x 2n matrices in the replica space, corresponding to the decoupling of the 5* pairs with the 
small transverse or small sum momentum, respectively. Then, in going from Eq. H2()|) to Eq. (|27|l , one integrates 
independently over fields d and c which leads to the appearance of the factor 1/2 in the latter equation. Some 
familiarity with the standard diagrammatics may help one to recognise that d corresponds to a diffuson - the 
sum of impurity ladder diagrams with a small transfer momentum, while c corresponds to a cooperon - the sum 
of ladder diagrams with a small sum momentum, Fig. 01 To make it a working analogy, one needs to use the 
fact that in Eq. H27|) Q is a slow varying field. To this end, one should expand the functional H27|) in gradients 
of this field. 
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B. Gradient Expansion 



The first step is to define a class of stationary, spatially-homogeneous (i.e. 'infinitely slow" ) fields for which 
the functional H27|) satisfies a saddle-point condition. By varying the functional with respect to the field Q, one 
finds 



ittvqQ = (r| Q\r) = Ir 



(30) 



This saddle-point condition is obviously satisfied at zero frequency uj and in the absence of the source field tj 
by any spatially- homogeneous matrix Q of a proper rank and symmetry that obey the condition = 1. This 
condition allows one to invert explicitly the matrix in the square brackets in Eq. (|30() : in the Fourier momentum 
space ^ = /2m — e-p becomes just a scalar variable of integration so that one finds 



g{p, c^=o) = 



' U + ^Q]. (31) 



After this, the corresponding integration over all momenta in the r.h.s. of Eq. (|30(l is straightforward in the pole 
approximation (justified provided that the week disorder condition, £f ^ci 3> 1, is fulfilled), when one replaces 
this integration by vq d^. 

The condition = 1 is convenient to parameterise as 

Q^U^AU, (32) 

where U belongs to a so called symplectic group Sp(2n). This means that U'^U = 1, and each ?7 is a "real 
quaternion" matrix, i.e. can be represented as 2n x 2n matrix with matrix elements of the structure a -I- 6i -I- 
cj -\- dk, where a, b, c, d are real n x n matrices and i, j, k are the so called quaternion units, represented as ir 
with T being the set of 2 x 2 Pauli's matrices in the "charge-conjugate" space defined by Eq. H15|l . Such a choice 
of U preserves the charge-conjugate symmetry of matrix Q, Eq. H29|) . However, each Q can be represented in 
the form (|32fl in infinitely many ways, as the "gauge transformation" U SU leaves Q unchanged for any S 

commuting with A = diag (1„, — 1„). Such a matrix 5* is a block-diagonal matrix S = diag ^S'„, S'„^ , with each 

Sn belonging to the symplectic group Sp (n), i.e. S* G Sp (n) (g) Sp (n). To make the representation (|32|) single- 
valued, one needs to factorise out matrices S (i.e. to fix the gauge). As a result of such a factorisation, one ends 
up with a coset space: U G Sp (2ri)/Sp (n) (g) Sp(n) for matrices U that makes a single- valued representation 
for Q in Eq. (|32|l (as matrices Q and U are now in one-to-one correspondence, Q belongs to the equivalent 
coset space, although its matrix representation would be obviously different). A general geometric structure of 
such a coset space is not difficult to understand and use (see, e.g., an excellent bool*^ addressed specifically to 
physicists). However, it won't be necessary since in what follows an explicit single- valued parameterisations of 
matrices U and Q will be used, which could be straightforwardly verified without any benefit of knowledge of 
differential geometry. 

Note that Q = A is a solution to the saddle point equation at nonzero uj. What is usually not stressed 
is that the saddle-point equation is satisfied not only by A but by any spatially-homogeneous diagonal matrix 
commuting with A with eigenvalues equal to ±1. This must be taken into account in considering solutions 
with replica-symmetry breaking, the latter being the necessary conditions of obtaining non-perturbative result 
within the replica approach22i2^. However, the existence of this wide class of the saddle-point solutions at 
Lu ^ is irrelevant for most applications and, in particular, for the expansion of the functional of Eq. H27() near 
the saddle point of Eq. I|3U|) . In such an expansion, the transverse modes with 7^ 1 can be neglected as, 
due to the first term in the functional H27|) they would incur a gap of order of 1/tciS S> 1 (for a homogeneous 
excitation) and would be thus suppressed. Then the first term gives just an irrelevant constant that could be 
neglected (it vanishes in the replica limit as Tr 1 cx n but such a term would be irrelevant even in a finite-n cr 
model.) Terms with a smaller gap of order to/ 6 would appear upon the expansion of the second, Tr In term. 
Keeping them is justified as long as loToi <C 1 which will be assumed throughout. The expansion of the Tr 
In term will also produce gapless gradient terms with the magnitude of order (l/T(.jS){qi)'^ in the momentum 
space: keeping them while neglecting the ^ 1 modes is justified as long as <C 1. Therefore, the two 
conditions of the applicability of the standard NLctM are 

WTel < 1 g£ < 1 , (33) 
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which should be taken together with the weak disorder condition, g ^ 1. Therefore, the standard NLcrM is 
appUcable for the description of the diffusive regime, at scales r ^ £ and times t ^ t^i. These conditions can 
be relaxed and the NLcrM can be generalisediSiiiiiiiSii& for considering also ballistic regime (r < £) but such a 
generalisation will not be considered in the present lectures. Now the expansion of the tr In term in Eq. (|27l) is 
straightforward by making use of the condition H32I) and applying the similarity transformation (omitting for 
now the source field r]): 



■ Trln 



- — Tr In 
2 



2toi 

- U 



Q 



LO ' 



A 



■ Tr In 



2t, 



-A 



-uku^ 



e>t]+|c/AL/t} 



(34) 



const - ^ Tr {GqU KU"^) + - Tr (g^U 



GoU 



where Go is Green's function (|3HI with Q = A and the symbol Tr implies the integration over r as well as the 
summation over matrix indices. The real part of the first term is an irrelevant n-proportional constant, while its 
imaginary part reduces to — (i[.i;Q;/4) Tr (Af/Af/^) = — (i[x)a/4) Tr (AQ), with the coefficient of proportionality 



a = AQmC/o(r,r) = A3m jGoip) 



' J df (^^ + l) = VTT . Thus, the first term in the expansion 



gives the gradientless term in the NLcrM functional Q ■ In the second term in the expansion H34|) , the operator 



^ = (— l/2rn)V^ — £f in the commutator 



^, W reduces to • V (n is a unit vector) so that Fourier- 
transforming it into the momentum space and performing the angular integration, one finds this term equal 
to 



{r\Go\r') (n-A)^, {r'\Go\r) {n ■ A\ 



d'^rd'^r' = 



4d ^ 



pA.q 



:^0,p+q*^q 



where A{r) = C/(r)VL/^(r) = —{'VU)W. Neglecting higher-order terms in q, one has Go,p+q 
tuting J2p ^ J d^i one reduces the above expression to 



2d 



E 



{^+iA)Aq ■ {^+iA)Aq 

ie + if 



^tr[Aq, Af = ^Tr(VQ)' 



8 



8 



Go,p; substi- 



(35) 



where D — WpTd/d is the diffusion coefficient, and the identity VQ = W [A, A] U was used, that follows from 
the representation H32|l and the definition of A. This completes the derivation of the NLcrM functional Ql, 
the matrix field Q obeying constraints Q or, equivalently, (|32|l . Including the symmetry-breaking factors, 
like a magnetic field, is straightforward as I will demonstrate below when deriving the NLcrM for interacting 
electrons. The source fields r) were left out of the expansion, but including them is simple: as seen from 
Eq. (|25|) . they enter on par with the term to A, so that their inclusion would result in additional gradientless 
terms, Tr rjQ. In a similar way, one can include any other source fields, e.g. those necessary for direct evaluation 
of the conductance^^. 

The expansion H34|l was made to the first non- vanishing powers in the small parameters 1)33(1 . The higher 
order terms, both in gradients and in the (symmetry-breaking) frequency may also be included. They are 
required in order to obtain tails of different distribution function or to describe the long-time relaxation in the 
framework of the replica NLcrMii but this goes far beyond present considerations. 

What is left to do, before turning to the interaction case, is to illustrate how the NLcrM works. The model is 
nonlinear because of the presence of the constraint ||2Jl - so such an illustration should show how the constraint 
may be explicitly resolved. From all the examples listed in the introduction, I will consider only one - that of 
mesoscopic fluctuationsii*^^i2^. 



V. MESOSCOPIC FLUCTUATIONS WITHIN THE NLaM 



The most famous of sample-to-sample (mesoscopic) fluctuations are, of course, the universal conductance 
fluctuation aS i "*"^ ! whose full distribution was first calculated within the replica NLcrMii. For the illustrative 
purposes of this section, it is more convenient to consider mesoscopic fluctuations of density of states (DoS), 
as this is the simplest quantity to calculate, whether diagrammatically-'^ or within the NLcrM^^. The main 
simplification within the NLcrM is that one does not need to introduce any source field, since all the moments of 



13 



the DoS fluctuations, e.g. variance, can be obtained by differentiating the partition function H28|) with respect 
to frequency lo: 



2 n^o yirn I duj^ 



(36) 



Here Sv = i/ ~ i/q, i^q is the disorder-independent average DoS, and 6q = 1/i'oL'^ is the mean level spacing. To 
get Eq. (|36|l . one first expresses v in terms of Green's functions. 



and using the representation of in terms of functional integrals, Eqs. H13 |l - H18(l one represents the above 
expression as 

^ ^ /l?^I?Ve'^^(^'")Tr^(r)At/;(r)- * ^^^^""^ 



27rL^ Zi J ^ ^ V / w ttU' du ' 

The expression for (^{Sv)^) is then obtained by the replication and the ensemble averaging of the product of 
two DoS, Eq. (|22l) . as described above. Then Z{uj) is finally expressed as the partition function for the NLcrM, 
Eq. H28|) . It is implied in Eq. H36|l that one considers the fluctuations in an open sample, where an inevitable 
infrared divergence is cut off at w '--^ i?T (which allows one to put w ^ 0) where the Thouless energy i?T = D/L"^ 
is inversely proportional to the diffusion time through the sample. For a closed sample, Eq. H3t)|l at finite uj can 
be used for calculating the correlation function {v[e + uj)v{e)) but for a; < <5o considerations should be more 
subtle (seeSSiS^) than what follows. 

The simplest task is to calculate Z[ijj) pcrturbatively. In order to make any calculation, one has to resolve the 
constraints (0), i.e. to express Q in terms of unconstrained matrices. For the purpose of constructing regular 
perturbation series, it is most convenient to parameterise Q as foUowsii: 

Q {t-W/2)K{t-W/2)-^ . (37) 

Here W is anti-Hermitian matrix, off-diagonal in the "retarded-advanced" space, labelled by + or — indices as 
described between Eqs. H21|l and H23() : 



^ - [W^^ j - {-B^ j (38) 

Here the 2n x 2n matrices B are unconstrained but for the fact that they obey the charge conjugation symmetry, 
Eq. 11511 . This means that they are "real quaternion" matrices, as described after Ea. (|32() : B — Wl2+'ib^Tf^, with 

being Pauli's matrices in the charge conjugate space and all [fi = 0, 1,2,3) being totally unconstrained 
n X n real matrices. The parametrisation H37|l makes Q automatically obeying both constraints in Eq. |(2Jl, 

= 1 and trQ = 0. The structure (|38|) makes W anticommuting with A, thus ensuring that Q is Hermitian. 
The functional integration in Eq. H28|) should now be performed over independent matrix components of B. A 
particular convenience of the parametrisation H37|l is in the fact that in the n = replica limit, the Jacobian of 
the change from Q{r) to B{r) is 1. To perform the integration, one should (again!) reduce it to Gaussian, by 
representing = J^o + {J^ — To)- Here the quadratic in W functional To is obtained by expanding the NLctM 
functional to the lowest order in W: 

= /tr(W)^d'^r = ^Tr|Vi?r (39) 

The w-proportional term is not included into Tq (although it should, if one were interested in frequency- 
dependent observables), since in Eq. H3f)|l one needs only quadratic in lo terms. They are obtained by expanding 
e"'-^"-^") in Eq. (|2Hl for Z{lo) in powers of {T — Tq) and further expanding Q = A{l + W +^W'^ + . . . according 
to Eq. I|38|) . and finally performing the functional integration 



/ VB{...)exp{~To) 
J VBexpi^Tf,) 



(40) 
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As this integration is Gaussian, the result would be a product of all pair averages. The later can be found 
immediately, as Eq. (|39|l ensures that the the functional average (|4U|) splits into a product of Gaussian integrals 
over independent matrix elements of B: 

{bU^Kmir'))^ - n(r - r')d'^''S,kS,m , (41) 

where the propagator n(r— r') obeys the diffusion equation, 7ri^o^V^n(r — r') — S{r — r'), with proper boundary 
conditions (open sample). For an L'^ cube, 11 can be defined via its Fourier transform as 

n(r-r') = — y Are'-J-^'-'-'). (42) 
TT Dq^ 

q 

Eq. H41() is relatively easy to derive but very inconvenient to use, since all matrices that enter the averaging 
should be explicitly split into their components in the retarded-advanced space, like in Eq. H38|) . and in the 
charge conjugate space, like in Eq. H29|) . But all these components can be assembled together once and forever 
by casting the averaging formula (|41|) into those (two) written explicitly in the full matrix space of Qii: 

tT{WrPWr' Fij^^ = n(r- r')|tr(APAi?+ - Pi?+) +tr(AP) tr(Ai?) - tr(P) tr(i?)| 

trfw^^ p\ trfWr' r)) = n(r - r') tr(APAi? - APAR+ - PR + PR+) (43) 



Here P and R are arbitrary matrices of the same rank as W (or Q). Equations (|43|l allow one to calculate any 
perturbative (or RG) expression. Let us use it for calculating Z{uj) up to and thus finding the DoS variance 

In the lowest perturbative order one needs to calculate the functional average ^ (•?\2)q, where J^^ is the 2"^^ 
term in Eq. JQ). Expanding Q, one finds that 



■Fi^ = — -, — Tr AQ = — — 



Tr (ivW^) + ^ Tr (ivW*) 



(44) 



since the zeroth term is an irrelevant n-proportional constant, and all the odd order terms vanish as tr = 0. 
It is convenient to put uj inside the Tr: in general, source fields that are necessary to calculate any other 
quantity but DoS are matrix field, and in a similar expansion they would have been inside the trace. Now to 
the lowest order 

Z{u:)^-^{Tr{u;W') TricuW'))^ . (45) 

One applies the averaging formulae H43|l by choosing all possible pairings and performing them one after another 
(in an arbitrary order), considering at each step all matrices W (but two active, i.e. participating in the chosen 
pairing) as fixed. There are only two different pairing schemes possible in Eq. H45|) : one either make two inter- 
trace, or two intra-trace pairings. In the first case, one applies first of the formulae H43() . with P = R = 1; then 
only the last term does not vanish and equals (tr 1)^ = 16n^. Making two such inter-trace averaging one finds 
the contribution proportional to which vanishes in Eq. (|35|l in the replica limit. Therefore, one needs only 
the intra-trace averages which is convenient to represent diagrammatically, Fig.[Sl(a). In such a representation, 
each trace corresponds to a (hatched) polygon with vertices representing either matrices W or the source field 
UJ (no more than one vertex in each polygon). At the first step of averaging, one applies the second of Eq. H43|l 
with P = ujW{r) and R = LLjW{r'). Taking into account a combinatorial factor of 2 (each W in the first trace 
can be paired with either of two in the second trace) and using that W is anti-Hermitian and commutes with 
A, one finds 



Combining this with Eq. H45|l and substituting the result into Eq. (|36(l . one recovers the well-known result^: 

((H^)_^MV^-^ (46) 
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a b c 

FIG. 5: One-loop (a) and two-loop {b and c) contributions to the DoS variance: hatched polygons with 2n vertices 
correspond to the expansion of tr(a;AQ) to the 2n"^ power in W , with external vertices representing the source field a;; 
hatched 2n polygons correspond to the appropriate expansion of tr(VQ)^; the wavy lines represent W pairings, i.e. the 
difi'uson propagators H42|l . in correspondence with the standard diagrammatic techniques, Fig.0] 

where g — is the dimensionless conductance of a cube of volume i'*, and Cd is a numerical 

coefficient dependent on the sample shape and dimensionality. 

Exactly the same technique could be applied for calculating the two-level correlation function, 
{v{£)v{£ + w)) /vq — 1. In this case u) plays a dual role of the frequency and the source field. The only 
technical change is that the TrujW'^ should be included into the zeroth approximation, Eq. (|39f) . which results 
in replacing Dq^ by Dq^ — iuj in the denominator of Eq. (|42|l . In the diffusive regime of relatively large uj, i.e. 
Lu > D/L^, the result H4()l) remains roughly valid, with substituting L^^ = {D/ujy^^ for L into the expression 
for g. However, this would fail in the two-dimensional case where for <^ L the coefficient C2 in Eq. (|21(l 
vanishes due to the analytic properties of the diffusion propagator n(a', q)^. In this case, the main perturbative 
contribution is due to the two-loop diagrams. Fig. [3(6, c). The hatched pentagon in diagram (6) corresponds to 
the 2"'^ term in the expansion (|44|1 . The hatched square in diagram (c) results from the average {T^Ty ^)^ with 

J-y''' obtained by expanding VQ^ to the 4"^ power of W. Such expansions give all the necessary coefficients 
corresponding to the Hikami boxes of the standard diagrammatic technique, so that one is left with calculating 
two-loop integrals containing 3 or 4 diffusion propagators. Note that in the diagrammatic technique the dia- 
grams b and c cancel each other while the only contribution results from the third possible 2-loop diagram^ 
describing the weak localisation correction to one of the diffusons. Such a diagram, on the contrary, vanishes 
within the parametrisation (|37|l while the identical result is given by the remaining two- loop diagrams (6, c). 

Finally note that the RG calculations are also based on the averaging formulae (|43|l which could also be 
appropriately modified to include cases of partially broken symmetry, either due to the presence of a magnetic 
field or magnetic impurities (unitary symmetry) or to the presence of spin-orbit scattering (symplectic symme- 
try). In all cases, the same averaging formulae are also applicable to the case when the NLcrM also encompasses 
interactions. 



VI. COULOMB AND PAIRING INTERACTIONS IN THE SIGMA MODEL 

The NLcrM works perfectly for perturbative problems and - in certain cases - does reasonably well even for 
non-perturbative problems within the model of non-interacting electrons in the presence of disorder. Since the 
inclusion of interactions beyond the mean-field approach is impossible within the supersymmetric method, the 
replica approach is still the most reliable tool for models with interactions. 

The original fermionic replica a model^i has been generalised by Finkelstein^ to include the interactions. 
The interest in this approach has been greatly enhanced by the recent discover}!^ of an apparent metal- 
insulator transition in 2D disordered systems in zero magnetic field. It is not at all clear whether the observed 
effects are, indeed, due the transitionS^iS, and if they are - whether such a transition is, indeed, driven by 
interactions. However, a possibility of having the metal-insulator transition in a 2D disordered interacting 
system (schematically described by the TOE model Q), while it is absent in a disordered noninteracting 
system where all states are localised, is intriguing by itself. Such a possibility is arguably a driving force in a 
considerable revival of interest in Finkelstein's cr model (see, e.g., Refs. 17 18. 19 20,.2L') . At the moment, there 
is no clear evidence whether Keldysh techniques employed in RefsiiSiiS would give any edge over the replica 
method described in these lectures. So all further considerations are also limited to the replica approach. 

Below I describe how to include the BCS pairing interaction and the Coulomb interaction (for simplicity, in 
a singlet channel only, as the inclusion of a triplet channel is technically almost identical) into the derivation 
of NLcrM given in section IIVI One inevitable complication is that such a derivation should be done in the 
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temperature technique, and here some leap of faith is necessary in order to come to any resuh within a 
relatively short paper. The resulting a model looks similar to that for the non-interacting model, although its 
usage is considerably mode difficult. Still, I give a relatively simple example of such a usage to derive some 
known result for the BCS model. Note that the derivation follows that in Refi^ which is similar in spirit to 
that in Ref.^^ but rather different from the original derivation in Refi. The derivation below is also made as 
parallel as possible to that of the NLcrM in the absence of interactions. 



A. Hubbard— Stratonovich decoupling 



Since the pair interaction, in contrast to the elastic scattering from impurities, does not conserve single- 
particle energy, the effective interaction functional should be dynamical. As usual, it is convenient to introduce 
imaginary time r implying everywhere the thermodynamic Gibbs averaging (together with the averaging over 
quenched disorder where applicable). Thus one considers thermodynamic Green's functions related (in the 
Matsubara frequency representation) to the retarded and advanced Green's functions of the previous sections 
by the standard procedure of analytical continuation^". Then the effective functional corresponding to the 
interaction term in Eq. ^ can be written as 



(47) 



where x = {r, r}, and all the fermionic fields are anti-periodic in imaginary time r with period 1/T, s — ("f , |) 
is the spin index, Vxx' = 6{t — t')Vc(t — r'), and Vc represents the Coulomb interaction. The fields are 
still assumed to have the same matrix structure, Eq. H15|l . as in the nonintcracting case, necessary to take 
into account both the slow modes (diffuson and cooperon) arising from the impurity scattering. After the 
replication necessary for treating the random potential, all the fields in Eq. (I47II naturally have the same replica 
index and the action would be just a sum of all the actions (|47|) with different replica indices. Similarly, the 
BCS functional describing the attraction in the Cooper channel can be written as 



Ao dxip'f{x)'ipi{x)ipi{x)ip'j{x), 



(48) 



where Aq is the BCS coupling constant. The total action is now the sum of 

The Hubbard-Stratonovich decoupling of the functional (|47|l - (|48|l is similar to that for the disorder functional 
(|26|l . i.e. some new matrix fields are introduced to decouple the quartic terms via the Gaussian integration like 
in Eq. H26|l . To perform this simultaneously with averaging over disorder, one needs to replicate all the fermionic 
fields as described after Eq. (|21|l . All the quartic interaction terms in the action are naturally diagonal in the 
replica indices. 

To allow for all the slow modes in the interaction functional ()47|l one should introduce three matrix fields: 
the fields $ and / to take account of a small-angle scattering (a singlet channel) and large-angle scattering, and 
yet another one that corresponds to the Coulomb repulsion in the Cooper channel. The last one would lead to 
the standard renormalization of the BCS attraction (see, e.g., Ref.^) and is not considered here; further, it is 
assumed that systems under considerations still have an effective attraction in the Cooper channel after such 
a renormalization. The interaction in the Cooper channel is then decoupled with the help of another matrix 
field, A. Then the decoupling has the following form: 



Jv<i>cxp l^-^Jdxdx' <^{x)V^J<^>{x') + i Jdx-iis{x)'^>{x)'i's{x) 
yp/expj-i Jdxdx' tr[/(x)V;;}/(x')] +^ Jdx^s{x)fss'{x)^s'{x) 
JvAexpl^~j J da;|A(a;)|^ -l-z Jdx [A(a;)l'|(a;)Wx(x) - A(a;)*|(a;)*|(a;) 



(49) 



After performing the Gaussian integration here, together with that in Eq. (|26|l . one obtains instead of the action 
(|27|l the following effective action that depends both on a (which is a new name for Q, the latter being reserved 
for a later usage) and on all the new fields: 



iclds 



— Tr CT Tr In 

8t 2 



^-i[—a + ^ + A + e]-f 



(50) 
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Here ^Fgoids is the quadratic in $, / or A part of the action H49|l : the operator e equals ir^dr in the imagi- 
nary time representation with T3 being the appropriate PauH matrix in the charge-conjugate space defined by 
Eq. H15|l (it becomes a diagonal matrix of fermionic Matsubara frequencies, e„ = 7r(2n -I- 1)T, in the frequency 
representation); Tr refers to a summation over all the matrix indices and to an integration over r and r (or 
summation over Matsubara frequencies in the frequency representation). The triplet channel, described by 
the field /, is quite important: in particular, it can lead to the delocalization in the presence of disorderiii. 
Nevertheless,, such effects will not be considered here and this term will be ignored from now on - mainly for 
simplicity (also, it does not contribute when considering effects of the Coulomb interaction in superconducting 
systems). The field a has the same structure as Q in Eqs. H29() and (|32(l . apart from explicitly including the 
2x2 spin sector and replacing the 2x2 retarded-advanced sector by the dependence on the imaginary time 
T (since it is diagonal in r, it becomes a matrix field in the Matsubara frequencies). In the absence of the 
fields $ and A, the saddle-point solution for the action H50|l is formally the same as in the zero-temperature 
techniques of section Hvl Eq. H32|l . but the matrix A has a non-unit structure in the Matsubara (instead of 
advanced-retarded) sector: 



cr = U^AU , A = diag {sgn e} . 



(51) 



The field $, corresponding to the singlet part of the Coulomb interaction, is diagonal in all the sectors. The 
"order-parameter" field A is Hermitian and self-charge-conjugate, diagonal in the replica indices and coordinates 
r and r, and has the following structure in the spin and time-reversal space: 



fX(2:)T3 



(52) 



where |A| and x ^tre the amplitude and the phase of the pairing field A(r,r), fi and t^^ are Pauli matrices 
{i = 0, 1, 2, 3 with fg = 1) that span the charge- conjugate and spin sectors, respectively. 

When $ and A are included, the saddle-point equation for the functional (|50|1 can be formally written 
similarly to Eq. H30() for the non-interacting zero-temperature case: 



-iTriya{r) 



(53) 



Now CT = A does no longer represent the saddle point solution for e ^ 0. Still, one can derive an effective 
functional by expanding the above Tr In within the manifold (|51f) in the symmetry-breaking field e -f A -f $ 
and in gradients of a, as has been done in the original works by Finkelsteini. An alternative is to make first 
a similarity transformation around A within the manifold (|51|l to find the saddle point solution in the 
presence of the fields. This can be formally done with the help of matrix Uq that diagonalizes the Hermitian 
field e + A + $: 



£ + A + $ = U^qXUo: 



Jsp = U^AUo . 



(54) 



Here Uq (which is yet undetermined but fixed) belongs to the same symmetry group that defines the manifold 
(|32|l in the absence of the interaction fields. By substituting Eq. H54|) into the saddle-point equation H53|) . one 
can easily verify that this is, indeed, a spatially- homogeneous solution, provided that Atci3> 1 which will always 
be the case for a dirty superconductor (Axei^l) or a dirty metal (dimensionless conductance g^l). 

Now one can perform the expansion of Tr In, Eq. (|50|) . in gradients of Q and in the symmetry-breaking fields 
represented by the eigenvalues A and the matrix Uq, Eq. (|54|l . by employing the following parametrisation: 



U^'AU. 



(55) 



where Q represents the saddle-point manifold in the metallic phase and a is obtained from Q by the same 
rotation (|54|) as agp is obtained from the metallic saddle point A. Therefore, Q is defined, as in the metallic 
phase, on the coset space S(2n)/S(n) (g) S (n) where, depending on the symmetry, S represents the unitary, 
orthogonal or symplectic group. In the noninteracting case of section liVI only the case of the symplectic 
group, corresponding (unfortunately) to the orthogonal symmetry, was described, but a generalisation to the 
other two cases is straightforward and described in almost any paper on the NLaM for disordered systems. 
The parameterisation H55() simplifies considerably all the subsequent derivations and leads to a new variant 
of the nonlinear a model for interacting systemsSl which can be more convenient for many applications than 
the original one^. After substituting cr — U^U^ AU Uq, Eq. H55|) . into Eq. (|50|l . one obtains the following 
representation for the Tr In term: 



<5^ = -iTrln{Go-i + U^li U^^] - i{U\U^)} , 
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where one can also include an external magnetic field with the vector potential A: 

^ 1 

£. = 7;-{p~eA)^-eF. 
Zm 



ZTc\ 



The expansion to the lowest powers of gradients and A is now straightforward and similar to that for the 
noninteracting case. It results after some calculations in the following effective action: 



^=^Tr[$V;->]+^5: |dr|A.P + ^Tr 
The long derivative in Eq. H56|l is defined as 



^ {dQf - \Q 



dQ = VQ+ Ao-ieAh,Q =doQ+[Ao,Q] 



where the matrix Aq is given by 



(56) 



(57) 



(58) 



and 9o = V — [zeAfs, . . .] is the long derivative (|57|l in the absence of the pairing field A. Both t/g and A 
should be found from the diagonalization of e + A + $, Eq. I|54|) . Although such a diagonalization cannot be 
done in general, it will be straightforward in many important limiting cases. For A = $ = 0, the field Aq 
vanishes, d ^ do, and A ^ e, so that the functional (|56|l goes over to that of the standard nonlinear cr model 
for non-interacting electrons. 

The action H56(l is most general in the current context. It allows one to develop a fully self-consistent approach 
to superconductivity of dirty metals in the presence of Coulomb interaction. However, any application requires 
a set of further simplifications. As a simple illustration, I will show below how to use the model in a dirty 
superconductor near the metal-superconductor transition in the absence of Coulomb interaction. 



B. Ginzburg-Landau Functional 

In the vicinity of the metal-superconductor transition one can expand the action H56|) (in the absence of 
the Coulomb field <&) in the pairing field A. A further simplification is possible in the weak disorder limit, 
g ^ 1, when the field Q can be integrated out to obtain an effective action for the A-field only. In the 
quadratic in A approximation, the kernel of this action will give, with due account for the disorder, an effective 
matrix propagator of the pairing field which governs properties of a disordered superconducting sample near 
the transition. ^ 

To integrate over Q, one splits the action (|56|l into T = J- + where T is the standard NLcrM functional 
as in the metallic phase. Then one makes a cumulant expansion, i.e. first expands e"^'^"'"'^^^ in powers of J- a, 
then performs the functional averaging with e~-^° (denoted below by (. . .)g) and finally re-exponentiates the 
results. Apart from the last step, this is exactly as described in section for the noninteracting NLcrM. The 
expansion involves only the first and second order cumulants since the higher order cumulants generate terms 
of higher order in A. Then the only contributions to the action quadratic in A are 



T^A] - /dr|A^p_!^(Tr(A-6)Q 



Expanding A and Ag to the lowest power in A and performing the functional averaging one finds^J* the action 
quadratic in A as follows: 

^efr[A] = ^E /dr-A:(r)(r|/C.|r')A^(r'), (60) 

UJ 

with the operator fC^j given by 



e(w-e)<0 I 
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Here Il'^^^{r,r') = (r|n'^'''|r') are the cooperon and difFuson propagators. At 
that of the NLcrM for non-interacting electrons, Eq. (|42(l . In general, 



the later coincides with 



(62) 



where C = —D {V — 2ieA) defines the cooperon modes; H'* is obtained from H'^ by putting the external vector 
potential A = 0. In the last term in Eq. (|61() . n|*^| (0) ^nj^^i (r, r); this term may be obtained by expanding in 

the cooperon propagator with the renormalised diffusion coefficient, 



C 



1 



|(0) 



which is a weak localisation correction to free cooperon propagator Ilj^ 
The summation over Matsubara frequencies in Eq. (|61|l yields 



(r,r'). 



^ , r /I \u;\-C^ 

K-uj = In h V - + — 

To M2 AttT . 



2 / 47rr 



(63) 



where To = ^co {B = 0) is the transition temperature of the clean superconductor in the absence of a magnetic 
field and "0 is the digamma function. The weak localisation correction is proportional to the coefficient given 

by 



AirT 



nvV ^ -DgH V2 ^ 47rr 







AttT 



AttT 



For Lo — Q the coefficient ao = au)={){T) can be simplified in the two limits: 



^'(1/2) 



E 



oo = < 



^(1/2) 



L:S> Lt, 



L<^Lt: 



(64) 



where Lt = y/D/T is the thermal smearing length. 

The instability of the normal state (i.e. a transition into the superconducting state) occurs when the lowest 
eigenvalue of the operator IC^ becomes negative. The eigenfunctions of this operator coincide with the eigen- 
functions of the cooperon operator C. The lowest eigenvalue of C is known to be Cq = DB/cjjQ, where (jjQ is the 
fiux quanta. This ground state cooperon eigenfunction corresponds to the lowest eigenvalue /Co of the operator 
ICui- The condition ICq — implicitly defines the line Tc{B) in the (T, i?)-plane where the transition occurs: 



To ^ \2 



Co 



aoCo 
4^' 



(65) 



The term in the r.h.s. of Eq. (|65|l describes a l/y-correction to the main result. This weak localisation is 
linear in the magnetic field B and vanishes as B ^ as expected (the Anderson theorem^^). In a nonzero 
magnetic field the weak localisation correction to Be is positive which has a very simple explanation. The 
superconductivity is destroyed by the magnetic field when <E'(^) > $0: where $(^) is the flux over the area with 
linear size of the order of the coherence length ^ ~ ^ D jT and $o is the flux quanta. The weak localisation 
corrections reduce D and thus ^. Therefore, one needs a stronger field to destroy the quantum coherence. The 
same reasoning explains the growth of Tc in a fixed magnetic field. 

Note finally that it would be straightforward to include the leading weak-localization corrections in all orders 
of In by calculating the Q - averages in Eq. H59|l via the renormalization group. This would lead to the 
renormalization of D in the cooperon propagator H62|l . thus changing the shape of the Tc{B) curve. However, 
the value of Tc(0) will again remain unaffected, since the superconducting instability is defined by the onset of 
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the homogeneous zero mode in the operator JC, Eq. H63|l . which does not depend on the value of the diffusion 
coefficient in the cooperon propagator. 

It is worth stressing that the Anderson theorem reflects certain properties of a model rather than those of real 
superconductors. If one allows for Coulomb interaction, then the critical temperature of the superconducting 
transition is no longer independent of disorder. Combined effects of the interaction and disorder lead to 
corrections to the transition temperature proportional (in a slightly simplified way) to g-^ \t? {TcT^x)^^'^^. 
For a relatively weak disorder, the system still remains superconducting at T ^ while for a sufficiently 
strong disorder the above corrections would suppress the superconducting pairing at any temperature and 
make a 2D system insulating at T — > 0. Such a mechanism of suppressing A by disorder gives a possible 
scenario for a widely observed superconductor-insulator (SI) transition in two-dimensional structures [58-61] 
which is most adequately described within the NLcrMLiS similar to those developed in this section. However, a 
widely accepted scenario for such a transition is based on very different, both in origin and in implementation, 
modelsSSiS^iSSiSSiSi where |A| remains finite and the superconductivity is suppressed by the loss of the phase 
coherence due to the fluctuations of the phase x of the order parameter. These "phase-only" models can also 
be derived", under certain parametrically controlled assumptions from the NLctM developed in this section. 
However, the scope of these lectures does not allow me to go into any further demonstrations. 

In conclusion, I have demonstrated in detail how to derive from scratch the NLaM for non-interacting 
electrons in disordered media, have shown some example of its perturbative usage, and in a much more succinct 
mode, have illustrated how to derive and use it in the presence of the pairing and Coulomb interactions between 
electrons. 
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